Metastable state in a shape-anisotropic single-domain nanomagnet 
subjected to spin-transfer-torque 
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We predict the existence of a new metastable magnetization state in a single-domain nanomagnet 
with uniaxial shape anisotropy. It emerges when a spin-polarized current, delivering a spin-transfer- 
torque, is injected into the nanomagnet. It can trap the magnetization vector and prevent spin- 
transfer-torque from switching the magnetization from one stable state along the easy axis to the 
other. Above a certain threshold current, the metastable state no longer appears. This has important 
technological consequences for spin-transfer-torque based magnetic memory and logic systems. 



Spin-transfer-torque (STT) is an electric current- 
induced magnetization switching mechanism that is 
widely used to switch the magnetization of a nanomag- 
net with uniaxial shape anisotropy from one stable state 
to the other— «. A spin-polarized current is injected into 
the magnet to deliver a torque on the magnetization vec- 
tor and make it switch. This has now become the staple 
of nonvolatile magnetic random access memory (STT- 
RAM) technology^. 

In this Communication, we show analytically that the 
spin polarized current can spawn a metastable state in 
the magnet, which can trap the magnetization vector and 
prevent it from switching. This happens only if the spin- 
polarized current is smaller than a certain value. Thus, a 
minimum current - which may be larger than the critical 
switching current - may be needed for fail-safe switching. 

Consider a single-domain nanomagnet shaped like an 
elliptical cylinder with elliptical cross section in the y- 
z plane (see Fig. [T]). The major (easy) and the minor 
(in-plane hard) axes of the ellipse are aligned along the 
z-direction and y-direction, respectively. Let 9(t) be the 
polar angle and <fi(t) the azimuthal angle of the magneti- 
zation vector in spherical coordinates. 

At any instant of time t, the energy of the unper- 
turbed nanomagnet is the uniaxial shape anisotropy en- 
ergy which can be expressed as^: 



E(t) = B(t)sin 2 9(t) + constant term, (1) 



where 



B(t) = B{<f>(t)) = ^ M 2 s n [N d _ xx cos 2 (P{t) 

+N d - yy sin 2 4>{t) - N d ^ zz ] . (2) 

Here M s is the saturation magnetization, Nd- xx , Nd- yy 
and Nd-zz are the x-, y- and z-components of demagne- 
tization factor—, and Q is the nanomagnet's volume. 

The magnetization M(t) of the single-domain nano- 
magnet has a constant magnitude but a variable orien- 
tation, so that we can represent it by the vector of unit 
norm n m (t) = M(t)/|M| = e r where e r is the unit vector 
in the radial direction. The other two unit vectors are 
denoted by eg and for 9 and <j> rotations, respectively. 
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FIG. 1. A nanomagnet shaped like an elliptical cylinder. The 
nanomagnet cross-sections are on the y-z plane. The magneti- 
zation direction can be rotated with a spin polarized current. 



The torque acting on the magnetization within unit 
volume due to shape anisotropy 

T E (i) - -§r x VE{6{t),4>{t)) 

= -{2B(t)sin9(t)cos9(t)}e c/> - {B Qe (t) sin9(t)}e e , 

(3) 

where 

B 0e {t) = B 0e m)) = Y M 2 n(N d - xx -N d _ yy )sin(2(b(t)). 

(4) 

Passage of a constant spin-polarized current I perpen- 
dicular to the plane of the nanomagnet generates a spin- 
transfer-torque that is given byi 

TsttW = s [c s (V) sin6(t) e e - b s (V) sin9(t) e ] , (5) 

where s = (hj1e)nl is the spin angular momentum 
deposition per unit time and r\ is the degree of spin- 
polarization in the current /. The coefficients b s (V) and 
c s (V) are volt age- dependent dimensionless terms that 
arise when the nanomagnet is coupled with an insulating 
layer as in an MTJ 6 and the spin-polarized current tun- 
nels through this layer. We will use constant values of 
b s (V) and c s (V) for simplicity^. Furthermore, we will as- 
sume b s (V) — 0.3 |c s (V)| and |c a (V)| = 1 to be in approx- 
imate agreement with the experimental results presented 
in Refs. d,©]. For 9 = 180° to 0° switching, c s (V) = +1, 
and for 9 = 0° to 180° switching, c^V) = -1, while 
b s {V) = +0.3 for both cases. 

The magnetization dynamics of the single-domain 
nanomagnet under the action of various torques is de- 
scribed by the Landau-Lifshitz-Gilbert (LLC) equation 
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-^T effW (6) 



where T e ff(t) — Te(*) + Tstt(*)j a is the dimension- 
less phenomenological Gilbert damping constant, 7 — 
is the gyromagnetic ratio for electrons, and 
My = fioM s Q.. Using spherical coordinate system, with 
constant magnitude of magnetization, we get the follow- 
ing coupled equations for the dynamics of 9(t) and <p(t): 



1 ^ a2 )^ir = T^ti- s ( c s(V)+ab s (V)) 



+ B 0e (t)}sin6(t) - 2aB(t)sin6{t)cos6(t)} (7) 



(l + « 2 ) 



2^ _ M 

di My 



{s(b s (V) ~ ac s (V))+aB Q e(t)} 
+ 2B(t)cos9(t)] (sin9 ^ 0). (8) 



Note that when the magnetization vector is aligned 
along the easy axis (i.e. 9 — 0°,180°), the torque due 
to shape anisotropy, Te(£) and the torque due to spin- 
transfer-torque, Tstt(£) both vanish (see Equations ([3]) 
and (0)), which makes d6(t)/dt as well as d<p(t)/dt equal 
to zero. Hence the two mutually anti-parallel orientations 
along the easy axis become "stable". If the magnetiza- 
tion is in either of these states, no amount of switching 
current / can budge it (because the spin-transfer-torque 
vanishes), which is why these two orientations are also 
referred to as "stagnation points" . Fortunately, thermal 
fluctuations can dislodge the magnetization from a stag- 
nation point and enable switching. We will now show 
that there can be a third set of values (9 3 ,<p 3 ) for 9(t) 
and (f>(t) for which both d0(t)/dt and d<f>(t)/dt will van- 
ish. 

We determine the values of (#3,^3) as follows. From 
Equations ([7]) and ©, by making both d6(t)/dt and 
d<fi(t)/dt equal to zero, we get 

2aB(fo)cos0 3 = B 0e (fo) - s c s (V) - as b s (V) (9) 
2B(4> 3 )cos9 3 = -aB 0e (<j) 3 ) + asc s (V) - sb s (V). (10) 

From the above two equations, we get Bo e ((j) 3 ) = sc s (V). 
If we put Bo e ((j) 3 ) — sc s (V) in Equation (J9j) or in Equa- 
tion (fTU|) . we get 2B((j) 3 )cos8 3 = —sb s (V). Accordingly, 
we can determine the values of (#3, 4> 3 ) as 



1 



(h/2e)r)Ic s (V) 
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(11) 



(12) 



Note that <p 3 depends on c s (V) while 9 3 depends on both 
c s (V) and b s (V). Neither depends on the Gilbert damp- 
ing factor a. 



In order to understand the physical origin of the state 
(9 3 ,4> 3 ), consider the fact that the total torque T e ff (t) = 
M(i) x H e ff(£) can be deduced from Equations and 
© as 

Teff (t) = {-2B(t)cos8(t) - sb s {V)} sin0(t) e 

+ {-B 0e (t) + sc s (V)}sin9(t)e e . (13) 

We immediately see that T e ff (t) vanishes when 6(t) = 9 3 
and <p(t) — cf> 3 . Hence there is no torque acting on the 
magnetization vector if it reaches the state 6(t) = 9 3 and 
4>(t) = <f) 3 at the same instant of time t. Thereafter, it 
cannot rotate any further since the torque has vanished. 
Unlike in the case of the other two stable states where 
both shape-anisotropy torque and spin-transfer-torque 
individually vanish, here neither vanishes, but they are 
equal and opposite so that they cancel to make the net 
torque zero. If the magnetization ends up in this orien- 
tation, then it will be stuck and not rotate further unless 
we change the switching current / to change the spin- 
transfer-torque. Since changing / can dislodge the mag- 
netization from this state, it is not a stagnation point 
unlike 9 — 0°,180°. Hence, we call it a "metastable" 
state. 

If b s (V) = 0, then 9 3 = 90°, which means that 
the magnetization will be stuck somewhere in the x-y 
plane perpendicular to the easy axis, if it lands in the 
metastable state. This plane is defined by the in-plane 
and out-of-plane hard axes. Note that when c s (V) is neg- 
ative, 03 is in the range —90° < cf> 3 < 0°, but when c s (V) 
is positive, 0° < <j) 3 < 90°. The quantity b s (V) cannot 
be negative^. Note also that when / = so that there is 
no spin-transfer- torque, 9 3 = 90° and 4> 3 can be any of 
the following values: {0°, 90°, 180°, 270°}. Consequently, 
the magnetization vector is either along the in-plane hard 
axis (y-axis) or the out-of-plane hard axis (x-axis) . Since 
these are obviously metastable states in an unperturbed 
shape-anisotropic nanomagnet, we call the state (#3,03) 
a "metastable" state. 

For numerical simulations, we consider a nanomagnet 
of elliptical cross-section made of CoFeB alloy which has 
saturation magnetization M s = 8 x 10 5 A/niAS, and a 
Gilbert damping factor a — 0.01. We assume the lengths 
of major axis (a), minor axis (b), and thickness (I) to be 
150 nm, 100 nm, and 2 nm, respectively. These dimen- 
sions (a, b, and I) ensure that the nanomagnet will consist 
of a single ferromagnetic domain 11,12 . The combination 
of the parameters a, 6, I, and M s makes the in-plane 
shape anisotropy energy barrier height ^32 kT at room 
temperature. With the dimensions (a, 6, and I) chosen, 
the demagnetization factors (Nd-xxiNd-y Vl Nd- Z z) turn 
out to be (0.947,0.034,0.019)ii. The spin polarization of 
the switching current is always assumed to be 80%. 

We assume that the magnetization is initially along 
the +z-axis, which is a stagnation point. Hence, at 
K, no switching can occur. However, at a finite temper- 
ature, thermal fluctuations will dislodge the magnetiza- 
tion from the stagnation point and enable switching. At 
room temperature, the thermal fluctuations will deflect 
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FIG. 2. Switching delays for a range of switching current 700 fiA - 40 mA producing spin-transfer-torque for c s (V) = — 1 
and different values of b s (V). The switching current is varied in steps of 10 jj,A. (a) b a (V) = 0.3. No switching occurs for the 
following switching current ranges: 2-2.05 mA, 2.33-2.49 mA, 2.83-3.26 mA, 3.69-5.09 mA, 5.6-10.24 mA, 11.41-24.51 mA. 
The switching failure is either due to the metastable state or the magnetization vector ending up in a state where it stops 
rotating and begins to oscillate around a mean orientation. The latter is a different state not discussed here, (b) b s (V) = 0.05. 
Switching succeeds for the entire range of switching current. Metastable state appears for b s (V) > 0.05. (c) b s (V) = 0. 
Switching succeeds for the entire range of switching current. 



the magnetization vector by ^4.5° from the easy axis 
when averaged over time^, so that we will assume the 
initial value of the polar angle to be 9 ini t = 4.5°. We 
choose the initial azimuthal angle 4>i n it as +90° because 
it is the most likely value in the absence of spin transfer 
torque. Similar assumptions are made by others^. We 
then solve Equations ([7J and © simultaneously to find 
9(t) and <p(t) as a function of time. Once 9(t) reaches 
175.5°, regardless of the value of fat), we consider the 
switching to have completed. The time taken for this to 
happen is the switching delay. 

Fig. [2] shows the switching delays versus switching cur- 
rent for different values of b s (V). The switching delay is 
'infinity' in some current ranges when b s (V) = 0.03 be- 
cause switching failed [see Fig. 2(a)| . However, beyond 
the current of 24.51 mA, switching always occurs within 
a finite time, meaning that the magnetization never ends 
up in the metastable state. Simulation results show that 
if the value of b s is small enough (< 0.05), the metastable 
state does not show up [see Figs. 2(b) and 2(c)| . 

The important question is why switching fails only for 
certain ranges of the current /, i.e. why does the mag- 
netization vector land in the metastable state for cer- 
tain values of / and not others? The answer is that 
starting from some initial condition (9i n it,4'init), the 
angles 6{t) and fat) must reach the values #3 and fa 
at the same instant of time t. This may not happen 
for any arbitrary /. Hence, only certain ranges of / 
will spawn the metastable state. It is also clear from 
Equation (ITT1) that above a certain value of /, there 
will be no real solution for fa since the argument of 
the arcsin function will exceed unity. This value will 
be hhreshoid = [efj, M^n (N d ^ xx - N d - yy )] / [hrjc s (V)]. 
By maintaining the magnitude of the switching current 
above hhreshoid, we can ensure that the magnetization 
vector will never get stuck at the metastable state. For 
the nanomagnet considered, hhreshoid — 32.7 mA, but 
switching becomes feasible at even lower current of 24.52 
mA since in the range [24.52 mA, 32.7 mA], the coupled 



9 and (^-dynamics expressed by Equations ([7]) and flSJ do 
not allow 9(t) and 4>{t) to reach #3 and fa simultaneously 
starting from (9 iniu fa m t)- 

Another important question is whether thermal fluctu- 
ations can untrap the magnetization from this state. To 
probe this, we solved the stochastic LLG equation^ in the 
presence of a random thermal torque. Fig. [3] shows the 
magnetization dynamics for a switching current of 24.51 
mA at room temperature (300 K). We observe that the 
magnetization gets stuck at a metastable state (with #3 — 
97.58° and fa = 335.87°) ~50% of the time, which means 
that roughly one-half of the switching trajectories inter- 
sect the metastable state and terminate there. The values 
of 03, fa are also the angles predicted by Equations (fTTj) 
and (fl~2"]) . thereby confirming that the metastable state 
indeed has the origin described here. Increasing the tem- 
perature to 400 K helps by decreasing the probability 
that a switching trajectory will intersect the metastable 
stated. What is important however is that i/the magne- 
tization vector gets stuck in the metastable state and the 
current remains on, then no amount of thermal fluctua- 
tions can dislodge it. In other words, this state is stable 
against thermal perturbations. 

Unfortunately, an analytical stability analysis is pre- 
cluded by the complex coupled 9-<p dynamics. Therefore, 
we performed numerical stability analysis while spanning 
the parameter space as exhaustively as possible. In all 
cases, we found the state to be stable against thermal 
agitations. 

We also notice that oscillations precede settling into 
the metastable state. This is due to coupled 9- and 
(^-dynamics governing the rotation of the magnetization 
vector, which causes some ringing. The metastable state 
appears in the switching current range 11.41 mA to 24.51 
mA since within this range, 9(t) and fat) can reach #3 
and fa simultaneously starting with the initial conditions 
{Oiniufanit}- If the initial conditions are changed, the 
range can change as welli^. 

Finally, one issue that merits discussion is what hap- 
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FIG. 3. Room-temperature (300 K) magnetization dynamics when the switching current is 24.51 mA (c 3 (V) = —1, b s (V) = 0.3). 
(a) Dynamics of 0(t). (b) Dynamics of 4>(t). (c) The trajectory traced out by the tip of the magnetization vector in three- 
dimensional space. The magnetization gets stuck at a metastable state with #3 = 97.58° and 03 = 335.87°. This plot was 
obtained by the solution of the stochastic Landau-Lifshitz-Gilbert equation in the presence of a random thermal torque to 
simulate the effect of thermal fluctuations 4 . This is one specific run from 10,000 simulations performed in the presence of 
thermal fluctuations that shows that the latter cannot untrap the magnetization from this state at room temperature. This 
happens for all the 10,000 simulations if thermal fluctuations are brought into play after the metastable state is reached. This 
shows that the state is stable against room-temperature thermal perturbations. 
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FIG. 4. Time evolution of 9(t) when the switching current 
is 20 mA (c s (V) = -1, b s (V) = 0.3). Thermal fluctuations 
were ignored but we assumed Omit = 4.5° and 4>mu = 90°, 
respectively. The switching current is turned off at 2 ns after 
the magnetization vector gets stuck at the metastable state 
with 9 3 = 95.74° and 3 = 341.25°. The dynamics shows 
that the magnetization vector relaxes to the easy axis because 
of shape anisotropy, but the final orientation is the wrong 
orientation along the +z-axis rather than the desired final 
orientation along the -z-axis. Therefore, switching fails. 



pens if the spin polarized current is turned off after the 



magnetization gets stuck. In that case, the torque due to 
shape anisotropy will take over and drive the magnetiza- 
tion to the easy axis. One expects that if O3 < 90°, then 
switching will fail since the nearer easy axis is the un- 
desired orientation, whereas if O3 > 90°, then switching 
should succeed because the nearer easy axis is the de- 
sired orientation. Equation (|12p dictates that 63 > 90° 
since b s (V) is always positive. Unfortunately, these sim- 
ple expectations are belied by the complex dynamics of 
magnetization. The out-of-plane excursion of the magne- 
tization vector causes an additional torque that depends 
on 83, 4>3. The torque can oppose the torque due to shape 
anisotropy. As a result, even when 83 > 90°, switching 
can fail since the magnetization reaches the wrong orien- 
tation along the easy axis (see Fig. @|. 

In conclusion, we have shown the existence of a new 
metastable magnetization state in a single-domain nano- 
magnet with uniaxial shape anisotropy carrying a spin- 
polarized current. If the magnetization gets trapped in 
this state, switching will fail with non-zero probability 
even in the presence of thermal fluctuations. This has 
vital implications for STT-RAM technology. 
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In this supplementary material, we provide additional data and plots for magnetization dynamics 
with corresponding descriptions to elucidate the physics and behavior of the metastable state. It 



' is well known that there is a minimum (threshold) current needed to switch magnetization with 
G ' 

spin-transfer-torque that depends on various magnet parameters such as shape anisotropy energy 

O 

barrier, etc. [Ref. 2 in the main Communication]. The current ranges that we have considered 

(N , 

here are above this minimum current. The purpose is to show that merely exceeding the minimum 

in 

current does not guarantee switching because of the existence of the metastable state. Some excess 
current is needed to guarantee fail-safe switching. 

\d . 

In the main Communication, we showed that there is a range of current for which switching 
can fail for one specific initial condition Qi n it = 4.5° and (pi n u = 90°. In this supplement, we show 
that failure can occur for other initial conditions as well. Thus, the metastable state is a threat for 
H ■ almost any initial condition. 

Another issue we explore in depth is whether turning off the current once the magnetization gets 
stuck at the metastable will ensure that the shape anisotropy takes over and eventually switches 
the magnetization to the correct state, albeit with a long switching delay. We present simulations 
to show that this is not the case, although one might intuitively think otherwise. Even if the 
magnetization gets stuck at the metastable state after crossing the in-plane hard axis (#3 > 90°), 
turning off the current at that point will not necessarily make it relax towards = 180° and 
successfully complete switching. The strong coupling between the 9- and ^-dynamics can cause the 
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magnetization to rotate backwards (owing to the strong counter-clockwise torque produced by the 
out-of-plane excursion of the magnetization vector for certain values of (f>) and actually end up in 
the wrong state. These results show that the metastable state can be a serious spoiler in STT based 
switching. 
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Figure SI: Switching delay versus switching current in the range 700 - 40 mA. Thermal fluctu- 
ations were ignored in the magnetization dynamics, but the initial values of 9 and <ft were assumed 
to be 4.5° and 90°, respectively. Any lower switching current will increase the delay to more than 
5 ns and is hence not considered. We used the parameters c s (V) = +1, b s (V) = 0.3 in this plot. 
The positive sign of c s (V) signifies 9 = 180° to 0° switching. Unlike when c s = — 1 (9 = 0° to 
180° switching), the switching delay always remains finite showing that the magnetization vector 
never gets stuck at the metastable state within this range of switching current for anti-parallel to 
parallel switching. Changing 4>init did not make any difference. This did not happen because there 
is no metastable state in this case, but because 9(t) and 4>(t) cannot reach the values #3 and ^3 
simultaneously within this current range. Thus, even though the metastable state exists, it does 
not hinder switching since the magnetization vector never visits this state during its sojourn from 
one stable state along the easy axis to the other. In other words, the switching trajectory does not 
intersect the metastable state. 
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Figure S2: Switching delay versus switching current in the range 700 fiA - 40 mA. Thermal fluctua- 
tions were ignored in the magnetization dynamics, but the initial values of 6 and <f> were assumed to 
be 4.5° and 0°, respectively. We used the parameters c s (V) = —1 and b s (V) = 0.3 in this plot. This 
corresponds to parallel to anti-parallel switching. No switching takes place for the following ranges 
of switching currents: 2.49-5.18 mA, 7.1-31.08 mA. The switching dynamics was simulated for every 
10 //A of switching current in this interval. Note that the ranges of current where switching fails 
are sensitive to initial conditions since these ranges are different from those in Figure 2 of the main 
Communication, where the initial condition was 4>init = 90° and Qi n n = 4.5°. 
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Figure S3: Switching delay versus switching current in the range 770 - 40 mA. Thermal fluctu- 
ations were ignored in the magnetization dynamics, but the initial values of 6 and <f> were assumed 
to be 7° and 90°, respectively. We used the parameters c s (V) = — 1 and b s (V) = 0.3 in this plot. 
No switching takes place for the following ranges of switching currents: 2.08 mA, 2.5-2.66 mA, 
3.26-3.91 mA, 4.82-7.88 mA, 9.45-20.79 mA. The switching dynamics was simulated for every 10 
yuA of switching current in this interval. Note once again that the ranges of current where switching 
fails are sensitive to initial conditions since these ranges are different from those in Figure 2 of the 
main Communication, where the initial condition was 4>imt = 90° and 6i n u = 4.5°. 
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Figure S4: Switching delay versus switching current in the range 770 mA - 40 mA. Thermal 
fluctuations were ignored in the magnetization dynamics, but the initial values of 9 and <f> were 
assumed to be 3° and 90°. We used the parameters c s (V) = —I and b s (V) = 0.3 in this plot. No 
switching takes place for the following ranges of switching currents: 1.97-2.01 mA, 2.24-2.35 mA, 
2.63-2.89 mA, 3.22-3.87 mA, 4.21-6.15 mA, 6.49-12.67 mA, 13.62-28.2 mA. The switching dynamics 
was simulated for every 10 fiA of switching current in this interval. Once again note that the ranges 
of current where switching fails are sensitive to initial conditions since these ranges are different 
from those in Figure 2 of the main Communication, where the initial condition was 4>init = 90° an d 

&init = 4.5°. 
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Figure S5: Time evolution of 9{t) and <p(t) when the switching current is 24.51 mA (c a (V) = — 1, 
bs(V) = 0.3) and the temperature is 400 K. The magnetization vector gets stuck at the metastable 
state and cannot be unstuck by thermal fluctuations at 400 K. This plot was obtained by the 
solution of the stochastic Landau-Lifshitz-Gilbert equation in the presence of a random thermal 
torque to simulate the effect of thermal fluctuations. Occasional small perturbations around the 
metastable state show up in both the dynamics of 6(t) and 4>{t) due to thermal fluctuations; however, 
such small fluctuations are not sufficient to dislodge the magnetization vector from the metastable 
state. Simulation over a long time duration (500 ns) was carried out to confirm this. This is one 
specific run from 10,000 simulations in the presence of thermal fluctuations that shows that the 
latter cannot untrap the magnetization from this state at 400 K. This happens for all the 10,000 
simulations if thermal fluctuations are brought into play after the metastable state is reached. 
However, if we consider thermal fluctuations from the very beginning, magnetization does not even 
reach the metastable state for ~52% of the simulations and switching occurs successfully, while for 
the rest ~48% of the simulations, the magnetization vector gets trapped at the metastable state 
and thermal fluctuations cannot untrap it thereafter. The failure probability at room temperature 
(300 K) was ~50% instead of ~48%, which shows that increasing the temperature slightly decreases 
the probability of a switching trajectory intersecting the metastable state, but in the event such 
intersection occurs and the magnetization gets stuck, the elevated temperature still cannot untrap 
it. (a) Dynamics of 9(t). (b) Dynamics of (j)(t). 
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Figure S6: Time evolution of 9{t) and (j>(t) when the switching current is 24.51 mA (c s (V) = — 1, 
bs(V) = 0.3). Thermal fluctuations were ignored but we assumed Qinit = 4.5° and (frinit = 

90°, 

respectively. The switching current is turned off at 2 ns after the magnetization vector gets stuck at 
the metastable state with #3 = 97.58° and ^3 = 335.87°. The dynamics shows that magnetization 
ultimately switches to the correct desired orientation along the easy axis {9 « 180°) because the 
torque due to shape-anisotropy takes over once the current is switched off; however, the switching 
takes too long since the shape anisotropy torque is much weaker than the spin-transfer-torque, (a) 
Dynamics of 9(t). (b) Dynamics of 4>(t). This is one case whether turning off the switching current 
resulted in successful switching. The other case, discussed in the main Communication, resulted in 
switching failure. Thus, turning off the switching current once the magnetization gets stuck does 
not guarantee successful switching. 



